Persistence and dynamics in ANNNI chain 
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We investigate both the local and global persistence behaviour in ANNNI (axial next-nearest 
neighour Ising) model. We find that when the ratio re of the second neighbour interaction to the 
first neighbour interaction is less than 1, P(t), the probability of a spin to remain in its original 
state upto time t shows a stretched exponential decay. For re > 1, P(t) has a algebraic decay but 
the exponent is different from that of the nearest neighbour Ising model. The global persistence 
behaviour shows similar features. We also conduct some deeper investigations in the dynamics of 
the ANNNI model and conclude that it has a different dynamical behaviour compared to the nearest 
neighbour Ising model. 
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I. INTRODUCTION 

The tendency of a spin in an Ising system to remain 
in its original state following a quench to zero tempera- 
ture has been extensively studied over the last few years 
and is well-established as an example of the phenomenon 
called persistence in dynamical systems [1]. Quantita- 
tively, persistence is measured by the probability P(t) 
that a spin does not flip upto time t. P(t) shows a power 
law behaviour, i.e., P(t) ~ t~ e , where is a new expo- 
nent not related to any other static or dynamic exponent. 
This phenomenon has been observed and studied quite 
extensively in Ising models with nearest neighbour inter- 
action in different dimensions [1-3]. 

Apart from such 'local' persistence, one can also study 
the 'global' persistence [4] behaviour by measuring the 
probability Pc{t) that the order parameter does not 
change its sign till time t. At the critical temperature, the 
probability that the individual spins will not be flipped 
till time t has an exponential decay, while the global per- 
sistence shows an algebraic decay: Pc?(t) ~ t~ e ° . The 
critical temperature of the Ising model in one dimension 
being zero, the global persistence becomes a quantity of 
interest together with the local persistence at zero tem- 
perature. Exact values of the exponents for the nearest 
neighbour Ising chain are known to be 9 = 0.375 [2], and 
9 G = 0.25 [4]. 

In this paper, our objective is to investigate the ef- 
fect of frustration on (local and global) persistence, by 
computer simulation of an ANNNI (axial next-nearest 
neighour Ising) chain [5] at zero temperature. This model 
has the Hamiltonian 



H — — Y^ =1 {SiSi+\ — nSiSi+2), 



(1) 



where Si is the spin (±1) at i-th site and k (> 0) is the 
parameter which represents the amount of frustration. 
We choose this particular model for our study because 
it is perhaps the simplest classical model with tunable 
frustration and because this model as such shows very 



interesting static and dynamic behaviour [5-8]. As re- 
gards static behaviour, the ground state is ferromagnetic 
for k < 0.5, antiphase (++ - - type) for k > 0.5 and 
highly degenerate for k — 0.5. On the other hand, zero- 
temperature dynamics using single spin-flip does not lead 
to the ground state for < k < 1 but does lead to the 
ground state for k > 1 [7] . 

Our main observation is that for < K < 1 there is 
no persistence (i.e., no algebraic decay of P(t)) while for 
K > 1 there is persistence albeit with a persistence ex- 
ponent different from that of the unfrustrated nearest 
neighbour Ising model (i.e., k = 0). We also claim, using 
some novel approaches, that as regards domain dynam- 
ics, the ANNNI model belongs to a different dynamical 
universality class for k > 1. 

In section II, we describe the model and the studies on 
local persistence for n < 1 and K > 1. In section III, the 
dynamics of the domains for k > 1 is analysed in detail. 
The global persistence behaviour is presented in section 
IV. Since the behaviour at k = 1 is unique, we have dis- 
cussed it in a separate section (section V). The results 
are summarised and discussed in section VI. 



II. LOCAL PERSISTENCE IN ANNNI CHAIN 



The Model 

We take an ANNNI chain of L spins in one dimension 
with periodic boundary condition in a random configura- 
tion (infinite temperature) and quench it to zero temper- 
ature. Thus, our updating rule is that a spin is selected 
randomly from the system, it is flipped (not flipped) if its 
energy is positive (negative), and it is flipped with prob- 
ability 0.5 if its energy is zero. L is always chosen to be a 
multiple of 4 to ensure complete antiphase ordering. We 
have used L = 8000 - 12000 (unless otherwise mentioned) 
and averaged the results over 10 3 to 10 4 configurations 
as per necessity. 
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Case I : k < 1 

For local persistence, in the entire range < k < 1 
the decay of P(t) is not algebraic (Fig. 1), but rather a 
stretched exponential, 

P(t) ~ cxp(-at fi ) (2) 

with a = 1.06 and (3 = 0.45. The exponent (3 as well as 
the coefficient a is found to be independent of n. There 
is no special behaviour at n — 0.5, the multiphase point 
for the ground state. 
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FIG. 1. The variation of P(t) for k < 1 are shown for k 
= 0.2,0.4,0.6 against t 0A5 . The curves fit fairly well to a lin- 
ear form in the log-linear plot indicating that P(t) fits to a 
stretched exponential form P(t) ~ exp( — 1.06t°' 4j ). The local 
persistence seems to be independent of k. Here the simula- 
tions have been done for L — 12000. 

The dynamical processes occurring in the system are 
as follows. In the first few dynamical steps ( - the pre- 
cise number depends on the system size) domains of size 
1 are removed and one is left with domains of size > 2. 
No domain wall is annihilated or created henceforth but 
the dynamics continues indefinitely as all domains of size 
> 2 are unstable. As a result, the system attains a 
non-equilibrium steady state but the equilibrium state is 
never reached and in a finite time all the spins in the sys- 
tem are flipped. This justifies the faster than power-law 
decay of P(t) in the system. To illustrate our argument 
further, we observe from simulation studies that after the 
first few iterations (typically 10 iterations for a system of 
L = 16000 spins) the number of domain walls per spin 
(say M/L) attains a constant value of 0.2795. This con- 
stant value of M/L is quite close to the most probable 
value of M/L which can be theoretically estimated easily. 
Under the constraint that each domain is of size > 2, the 
number of configurations in a system of L spins and M 
domains, is 



and the maximum value of this quantity is for 
Case II : k > 1 

For local persistence, in the range k > 1, the behaviour 
of P(t) agrees well with an algebraic decay 

P(*)~t-*' (5) 

at large t (Fig. 2). However, 8' shows a weak dependence 
on the time interval over which it is calculated indicating 
that there is a correction to the scaling. We have calcu- 
lated the slopes in the log-log plot over different intervals 
of time and found that 8' increases slowly with time ini- 
tially but at large times (t > 10000) attains convergence 
to 8' = 0.69 ± 0.01. This value is considerably different 
from that of the ferrromagnetic Ising model. The value of 
8' remains the same for all k except at the point 1/k = 0, 
where the system breaks up into two independent sub- 
lattices with nearest-neighbour antifcrromagnctic inter- 
actions, and 8' becomes equal to 9. The dynamics here 
is such that the patches of antiphase state (e.g., ++ - 
- ) grow in size gradually and the equilibrium (lowest- 
energy) configuration of antiphase is reached eventually. 
Similar dynamics also prevails in the ferromagnetic state 
(k = 0) , as the regions of up (or down) spins grow in size 
and ultimately fill up the entire system. However, one 
must note that inspite of the similarity in dynamics, the 
persistence exponents are indeed different. In the next 
section we investigate the domain dynamics for n > 1 in 
greater detail. 
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FIG. 2. The variation of P(t) and Pa(t) are shown for k = 
2.0. Both have a power law decay, the best fit straightlines 
in the log-log plot for large t show that P(t) has an expo- 
nent ~ 0.69 while -Pg(£) has an exponent ~ 0.50 (see text for 
the details in calculating the exponent 8'). The curves are 
the same for all k > 1. Here the simulated system sizes are 
L = 10000 and L = 8000 for P(t) and P G (t) respectively. 
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III. DOMAIN DYNAMICS FOR n > 1 

In some earlier studies [6,7], it has been concluded that 
the ANNNI model belongs to a different dynamical uni- 
versality class compared to the nearest neighbour Ising 
model. However, precise values of exponents for other 
dynamical processes e.g., dynamical exponent, are not 
available for the ANNNI model. In the nearest neigh- 
bour ferromagnetic Ising case, the dynamical exponent 
z = 2. Here we have attempted to estimate the value of 
z', the dynamic exponent for the ANNNI chain. 

To estimate z' directly, we have calculated the lifetimes 
r defined as the time upto which the dynamics continues 
in the ANNNI model as a function of the system size L 
(Fig. 3). In the nearest-neighbour Ising model, this time 
varies as L z , while for the ANNNI model we find 



with z 1 ~ 2.15. 
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FIG. 3. Lifetime r (scaled down by 10 5 ) is shown against 
the system size L. r increases as L 215 . The value of k — 2.5 
here. 

The dynamics in the nearest-neighbour Ising model 
can be viewed as a random walk of the domain walls 
which separate regions of up-spins and down-spins. The 
domain walls annihilate each other when they meet and 
their number M reduces with time as i" 1 / 2 with the value 
of z equal to 2 as mentioned above. On the other hand, in 
the ANNNI model, the domain dynamics is different and 
the domain walls have a different role in the antiphase 
state. In order to compare the dynamics in these two sys- 
tems, it is therefore useful to regard the domain walls in 
the non-frustrated (k = 0) Ising model as "defects" and 
compare its dynamics with that of an analogous quan- 
tity in the ANNNI model. We have adopted a number 
of ways to estimate the defects in the ANNNI model (for 
k > 1) as a function of time to get an estimate of the do- 
main decay exponent. In the following we briefly describe 
these methods. 



First we note that the excitation energy AE(t) (devi- 
ation from the ground state energy) of any state in the 
nearest neighbour ferromagnetic Ising model is identical 
to the number of domains (apart from some multiplica- 
tive constant), and as the system relaxes to its equilib- 
rium ground state, AE also shows a decay AE <~ t~ 1//2 . 
In the ANNNI model also, we study the relaxation of en- 
ergy by computing AE which shows the behaviour (Fig. 
4) AE ~ t- 1 ' 7 -' with z' ~ 2.40. We claim that z' is the 
dynamic exponent in the ANNNI model. 
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FIG. 4. The time decay of M' (number of domains of length 
other than 2; curve maked I), of "defect" sites D (II) and of 
excitation energy AE (III) are shown against time t (k = 2.5). 
All the simulations have been done for systems with L = 1000 
spins. 

Secondly, a direct measure of the defects defined to be 

simply the spins that do not belong to +H type states 

has also been done (e.g., in the distribution +H h++ 

- -, the 7-th spin is a defect, while in +H h+++ the 7- 

th and 8-th spins are defects.) Our simulation shows that 
the number of such defects, D, also decreases as i -1 / z 
for all k > 1 with z 1 ~ 2.30. 

A third way to estimate z' is to count the number of 
domains which arc not of length 2, i.e., domains which 
do not satisfy the antiphase configuration. As men- 
tioned earlier, in [7] , the dynamics of domains of different 
lengths were studied from which it was concluded that 
the ANNNI chain belongs to a new universality class. 
Here we count all domains of size other than two. We 
find a power law variation of this quantity (say M') again 
with exponent 1/z' where z' ~ 2.30. All the above results 
are shown in Fig. 4. 

From all the above estimates, we conclude that the dy- 
namic exponent exponent is z' = 2.3±0.1 for the ANNNI 
model in one dimension. 

It is possible, in principle, to calculate the value of z' 
also from the relation 



P(t -» oo, i) ~ L 



-z'e' 



(6) 
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when the dynamics has stopped altogether [9]. How- 
ever, such an estimate would require simulations upto 
L z number of Monte Carlo steps for a system of size L. 
We could get these data only for systems with L < 1000 
which is insufficient to give a good scaling of P(t — > oo, L) 
with L. 

IV. GLOBAL PERSISTENCE 

We shall now present our results on global persistence. 
In this context one needs to be careful in defining the 
order parameter for different values of n. For k < 0.5, 
the global persistence may be calculated in terms of the 
magnetisation as for the ferromagnetic nearest- neighbour 
Ising model and the global persistence shows an exponen- 
tial decay : 

P G (t) ~ exp(-O.lOi). (7) 

In the entire range 0.5 < k < oo, one needs to define the 
antiphase order parameter. It is not easy to put forward 
such a definition in a straightforward manner. One may 
be tempted to define the sublattice magnetisation 

Si + S2 — S3 — S4 + 1S5 + — S7 — Sg + • • • 

as the order parameter but such a definition docs not 
work because the antiphase has a four-fold degeneracy. 
We propose instead to consider the magnetisations over 
four sublattices : 

m a = ^Lt 1 S a+i3 - a = 1,2,3,4 (8) 

as the order parameter. The global persistence function 
Po{t) is evaluated for each sublattice and an average over 
the four sublattices is taken. We have verified that for 
k = 0, with this definition of order parameter one gets a 
power-law decay of Pc{t) with the known exponent 0.25. 



For 0.5 < k < 1, Po(t) again has an exponential decay 
P G (t) ~ exp(-0.16i). (9) 

which is faster than that observed in the region < k < 
0.5 (Eq. (7)). The results for global persistence for n < 1 
are shown in Fig. 5. 

For k > 1 there is algebraic decay of the global persis- 
tence (Fig. 2) 

p G (t)~t- e '°, (10) 

with B' G ~ 0.5, which is different from the exponent for 
the nearest-neighbour Ising model in one dimension. It 
should be noted that for n — the global persistence 
behaviour could be analysed exactly [4] because the dy- 
namics is essentially that of some independent random 
walkers which annihilate each other. There is no such 
simplification here, as the j-th spin in the a-th sublat- 
tice, namely S a+ 4j is independent of the neighbouring 
spins in the same sublattice, but depends on the j-th 
spin in the other sublattices. 



V. DYNAMICS AT k = 1 



The point k = 1 signifies a dynamical transition point 
with different dynamical behaviour on its two sides and 
the behaviour here is different from that for lower or 
higher k values. The fraction of persistent spins P(t) 
decays as a stretched exponential (Fig. 6) 
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FIG. 5. The variation of Pa{t) for k < 1 are shown. The 
labels indicate the values of k. The exponential decay is 
accompanied by different exponents for < n < 0.5 and 
0.5 < k < 1.0. The system size is L = 12000 here. 



P(t) - cxp(-1.92t - 21 ) (11) 

but this decay is slower than a similar decay (Eq. (2)) for 
< k < 1. As regards the behaviour of domain decay, 
we have found that (i) the excitation energy AE, (ii) the 
number of defects D (the sites that do not belong to the 

pattern +H ), (iii) the number of domains (M') which 

are not of length 2, all decay as t~ 1 / z and the lifetime 
r also grows with L as L z with z" ~ 3. 

It is interesting to note that the persistence P{t) does 
not decay algebraically which one would expect from the 
fact that the dynamics does lead to the equilibrium state 
although at a rate much slower than that for n > 1 [7]. 
This is again, we believe, a special feature related entirely 
to the persistence phenomena. 
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FIG. 6. The variation of P(t) and P G (t) for k = 1 are 
shown. Both have a stretched exponential decay, the best 
fit curves shown in the figure are 4.39 exp(— 1.69i ' 25 ) for the 
global curve and 6.52 exp( — 1.92f°' 21 ) for the local curve. Data 
correspond to a L = 8000 system. 

The dynamical transition point k = 1 is marked by a 
stretched exponential decay of the global persistence also 

P G (t) ~ exp(-1.69i - 25 ), (12) 

as shown in Fig. 6. 



words, the behaviour changes sharply as soon as a slight 
amount of next nearest-neighbour interaction is added 
to nearest-neighbour interaction, or a slight amount of 
nearest-neighbour interaction is added to next nearest- 
neighbour interaction. 

It should be mentioned that three regimes of persis- 
tence (namely, exponential, stretched-exponential, and 
algebraic) have also been studied earlier, although in 
somewhat different contexts [10]. It should also be 
pointed out here that all the dynamical features obtained 
here are for the single spin- flip Glauber dynamics. Other 
types of dynamics, as considered in [6,7] could lead to 
different behaviour. 

Before we conclude we must also mention that for se- 
rial updating (where the sites are visited serially and up- 
dated) the results presented here remains qualitatively 
the same in the sense that the indices remain the same 
upto the accuracy of the simulation. Also, for n > 1, 
one obtains the same values of the exponents even if one 
starts with a ferromagnetic state instead of a random 
state. 
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VI. SUMMARY AND CONCLUSIONS 

In summary, we have studied the persistence behaviour 
for the ANNNI chain for all possible values of n, the para- 
mater governing frustration in the model. Three regions 
of different behaviour of local persistence are obtained 
as k is varied. The fraction of persistent sites shows a 
stretched exponential decay for k < 1. At k = 1, it is 
also stretched exponential with a different exponent. For 
k > 1, algebraic decay of the persistence probability is 
observed with exponents different from that of the un- 
frustrated case. An estimate of the dynamical exponent 
z' using several approaches has been made from which we 
conclude z 1 ~ 2.30. This confirms that the ANNNI model 
belongs to a different universality class from that of the 
unfrustrated nearest neighbour Ising model for which the 
corresponding exponent has a value 2.0. When global 
persistence is considered, we again observe different types 
of behaviour of the persistent probability in the above 
three regions. The behaviour of both local and global 
persistence is unique at k = 1, the dynamical transition 
point. One should note that there is thus some abrupt 
change in behaviour (both of local and of global persis- 
tence) at the points k = 0, k = 1 and 1/k = 0. In other 
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